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ABSTRACT 

A  quasi-three-dimensional  thin-layer  Navier- 
Stokes  analysis  was  used  to  predict  heat  transfer 
to  rough  surfaces.  Comparisons  are  made  between 
predicted  and  experimental  heat  transfer  for  tur¬ 
bine  blades  and  flat  plates  of  known  roughness. 
The  effect  of  surface  roughness  on  heat  transfer 
was  modeled  using  a  mixing  length  approach.  The 
effect  of  near-wall  grid  spacing  and  convergence 
criteria  on  the  accuracy  of  the  heat  transfer  pre¬ 
dictions  are  examined.  An  eddy  viscosity  mixing 
length  model  having  an  inner  and  outer  layer  wets 
used.  A  discussion  of  the  appropriate  model  for 
the  crossover  between  the  inner  and  outer  layers 
is  included.  The  analytic  results  are  compared 
with  experimental  data  for  both  flat  plates  and  tur¬ 
bine  blade  geometries.  Comparisons  between  pre¬ 
dicted  and  experimental  heat  transfer  showed  that 
a  modeling  roughness  effects  using  a  modified  mix¬ 
ing  length  approach  results  in  good  predictions  of 
the  trends  in  heat  transfer  due  to  roughness. 

N  omencla  t  ure 


Prt-  Turbulent  Prandlt  number 
p  -  Pitch  between  roughness  elements 
Re  -  Unit  Reynolds  number 

5  -  Reference  surface  area 
Sj  -  Total  frontal  area 

s  -  Surface  distance 

St  -  Stanton  number  based  on  inlet  conditions 

U  -  Inviscid  velocity 

ur  -  Friction  velocity 

y  -  Distance  normal  to  surface 

y+  -  Normalized  distance,  {yUe/v)y/C / /2 

6  -  Full  boundary  layer  thickness 

k  -  Von  Karman’s  constant 

A  -  Roughness  density  parameter 

p  -  Dynamic  viscosity 

pt  -  Turbulent  eddy  viscosity 

v  -  Kinematic  viscosity 

vt  -  Turbulent  eddy  viscosity,  {pt/p) 

p  -  Density 

u;  -  Vorticity 


A/  -  Frontal  area 
A3  -  Windward  wetted  area 
A"1"-  Damping  coefficient 
Cj  -  Friction  factor 

D  -  Diameter  of  hemispherical  roughness 
h  -  Roughness  height 
/ieq-  Equivalent  roughness  height 
L  -  Distance  between  roughness  elements 


l  -  Length  scale 


Subscripts 

c  -  Edge  of  boundary  layer 
m  -  Measured 
p  -  Predicted 
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INTRODUCTION 

The  effect  of  inaccuracies  in  predicting  turbine 
blade  temperatures  significantly  affects  turbine 
durability.  As  discussed  by  Graham*1),  a  50°C 
underprediction  of  blade  temperature  can  decrease 
service  life  by  a  factor  of  two.  Taylor*2)  measured 
the  surface  roughness  of  in-service  turbine  blades. 
His  results  showed  that  for  typical  Reynolds  num¬ 
bers  the  surfaces  are  not  hydraulically  smooth. 
Since  actual  turbine  blades  are  rough,  heat  transfer 
predictions  assuming  smooth  surfaces  can  signifi¬ 
cantly  underpredict  external  heat  transfer.  This 
in  turn  could  result  in  underprediction  of  the 
blade  surface  temperatures.  Both  two-dimensional 
boundary  layer  analyses,  such  as  the  STANS  code 
of  Crawford  and  Kays*3),  and  Navier-Stokes  anal¬ 
yses,  such  as  those  of  Chima*4)  and  Hah*5),  have 
been  used  to  predict  smooth  turbine  blade  heat 
transfer.  Models  to  account  for  surface  roughness 
in  boundary  layer  analyses  have  been  developed 
by  Cebeci  and  Chang*6),  and  by  Taylor,  Coleman, 
and  Hodge*7).  Navier-Stokes  analyses  have  the  ad¬ 
vantage  over  boundary  layer  analyses  in  that  they 
can  be  used  to  predict  blade  heat  transfer  beyond 
the  point  of  separation.  The  present  work  dis¬ 
cusses  the  results  obtained  when  a  surface  rough¬ 
ness  model  was  incorporated  into  the  quasi-three- 
dimensional  thin-layer  Navier-Stokes  analysis  de¬ 
veloped  by  Chima*4). 

There  are  two  different  philosophies  with  re¬ 
gard  to  modeling  surface  roughness.  The  approach 
taken  by  Taylor  et  al.  *7)  explicitly  accounts  for 
roughness  geometry.  This  approach  accounts  for 
blockage  due  to  the  roughness  elements,  the  in¬ 
creased  drag  due  to  their  presence  in  the  flow  field, 
as  well  as  a  correlation  to  account  for  additional 
heat  transfer  between  the  fluid  and  the  discrete 
roughness  elements.  Consequently,  this  approach 
requires  detailed  knowledge  of  the  surface  rough¬ 
ness,  as  well  as  empirical  correlations  for  the  drag 
and  heat  transfer  to  the  roughness  elements.  Hosni, 
Coleman,  and  Taylor*8)  have  used  this  model  in  a 
boundary  layer  analysis  to  predict  heat  transfer  for 
a  rough  flat  plate  over  a  wide  range  of  Reynolds 
numbers.  The  other  philosophy,  as  expressed  by 
Cebeci  and  Chang*6),  is  to  model  the  roughness 
effect  based  on  the  offset  in  the  near-wall  velocity 
profile.  In  this  model  the  mixing  length  is  modified 


to  account  for  the  effective  height  of  the  roughness. 
The  effective  height  is  a  function  of  both  roughness 
density  and  shape.  Computationally,  the  model  of 
Cebeci  and  Chang  is  significantly  easier  to  imple¬ 
ment  into  a  Navier-Stokes  analysis  than  the  more 
detailed  model  of  Taylor  et  al.  AJso,  Hosni  et  al.*8) 
reported  that  their  boundary  layer  calculations,  us¬ 
ing  the  model  of  Taylor  et  al.,  were  done  using  250 
points  within  the  boundary  layer.  This  is  far  more 
than  is  practiced  to  use  in  a  Navier-Stokes  analysis 
for  turbine  blade  heat  transfer.  Consequently,  the 
model  of  Cebeci  and  Chang  was  implemented  in 
the  Navier-Stokes  analysis. 

Navier-Stokes  heat  transfer  predictions  can  be 
affected  by  the  near-wall  grid  spacing.  The  sen¬ 
sitivity  of  the  heat  transfer  predictions  for  both 
smooth  and  rough  surfaces  will  be  examined.  Be¬ 
cause  roughness  promotes  very  early  transition, 
heat  transfer  to  rough  surfaces  can  be  assumed 
to  be  fully  turbulent.  The  effect  of  different  ap¬ 
proaches  to  turbulent  eddy  viscosity  modeling  in 
the  leading  edge  stagnation  region  is  discussed. 
Comparisons  are  made  with  experimental  data  to 
demonstrate  the  suitability  of  the  roughness  mod¬ 
els  for  predicting  turbine  blade  heat  transfer.  The 
data  of  Hosni,  Coleman,  and  Taylor*8)  for  flow  over 
a  roughened  flat  plate  was  used  for  comparisons. 
Comparisons  were  made  for  the  heat  transfer  on 
turbine  blades  using  both  the  data  of  Liebert*9) 
and  the  data  of  Blair  and  Anderson*10). 

METHOD  of  ANALYSIS 

The  roughness  model  of  Cebeci  and  Chang*6) 
was  incorporated  into  the  Navier-Stokes  analysis  of 
Chima*4)  by  modifying  a  version  of  the  Baldwin- 
Lomax  *n)  turbulent  eddy  viscosity  model.  In  this 
model  there  is  an  inner  and  outer  region.  The  tur¬ 
bulent  eddy  viscosity  in  the  inner  region  is  given 
by: 

(Mt)iNNER  =  Pi*  M  (!) 

The  length  scale,  £,  is  given  by: 

e  =  Ky{  l-e"«'+'A+)  (2) 

Where  y+  =  yv/xiT  and  A+  =  26.  In  the  outer 
region  the  turbulent  eddy  viscosity  is  given  by: 

(Mt)ouTER  =  C'outer/’^’wake^’kleb  (3) 
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Where  Couter  =  0.0269,  and  Fwake  is  given  by: 


fwAKE  =  ™*n(yMAxf’MAX,-25yMAXt>DIF/J'MAx) 

(4) 

The  term  uDjF  is  the  maximum  velocity  along 
the  near  normal  grid  line  out  from  the  surface.  The 
quantities  i/max  and  Fmax  are  found  from: 

m = »|w|(i  -  «-»+/A+)  (5) 

J/max  is  the  value  of  y  for  which  F( y)  is  a  maxi¬ 
mum.  .Fkleb  is  given  by: 

^kleb  =  [l  +  0.004(y/yMAx)6]  (0) 

For  rough  surfaces  Cebeci  and  Chang  use  an 
offset  distance,  Ay,  which  is  added  to  the  physi¬ 
cal  normal  distance.  This  was  implemented  by  re¬ 
placing  y  by  y  +  Ay  in  the  above  equations.  For 
example,  equation  2  for  a  rough  surface  becomes: 

£  =  *(y  +  Ay)  (l  -  e~^++^/A+)  (7) 

At  the  wall  y  =  0,  but  £  >  0  when  the  surface  is 
rough.  Therefore,  fit  >  0  at  the  wall.  The  heat 
transfer  predictions  were  made  using  Jin  effective 
wall  thermal  conductivity.  Prt  was  assumed  equal 
to  0.9  for  the  determination  of  the  effective  wall 
conductivity. 

Cebeci  and  Chang  give  Ay  as: 

Ay+  =  0.9(yfh+l  -  h+e-</9)  (8) 

Cebeci  and  Chang  state  that  equation  8  is  valid  in 
the  range  4.535  <  <  2(300.  Figure  1  shows  the 

variation  of  Ay+  with  h+q.  When  the  equivalent 
roughness  height,  h+qy  equals  4.535,  Ay+  is  zero. 
This  automatically  accounts  for  the  fact  that  small 
roughnesses,  (/i+  <  5),  are  hydraulically  smooth. 
For  large  values  of  the  slope  of  the  curve  pro¬ 
gressively  decreases.  Consequently,  the  effect  of  in¬ 
creasing  roughness  height  becomes  less  significant 
as  the  roughness  becomes  large. 

In  the  boundary  layer  analysis  of  smooth  walls 
the  length  scale,  £,  is  held  constant  when  /cy  ex¬ 
ceeds  0.0865.  In  the  Navier-Stokes  analysis  the 
boundary  layer  thickness,  5,  is  not  clearly  defined, 
and  is  not  used  in  the  Baldwin-Lomax  turbulence 


model.  For  rough  surfaces  the  turbulent  eddy  vis¬ 
cosity  in  the  outer  region  was  calculated  in  two 
ways.  First,  y  was  unmodified  in  the  outer-region 
eddy  viscosity  equations.  In  the  second  way,  y 
was  replaced  by  y  +  Ay  in  the  outer-region  eddy 
viscosity  equations.  The  consequence  of  these  as¬ 
sumptions  on  the  rough  surface  heat  transfer  will 
be  discussed. 

The  Baldwin-Lomax  turbulence  model  has  an 
abrupt  crossover  model.  In  this  model,  when 
(fit) outer  >  (^t)lNNER,  f^t  is  taken  as  the  outer 
value.  Granville^12)  discusses  two  other  crossover 
models  which  blend  the  inner  and  outer  regions. 
These  two  models  are: 
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The  effect  on  the  predicted  heat  transfer  using  each 
of  these  models  will  be  discussed. 

RESULTS  and  DISCUSSION 
Heat  transfer  predictions  are  known  to  be  sen¬ 
sitive  to  the  near-wall  grid  spacing  in  a  Navier- 
Stokes  analysis.  This  sensitivity  could  be  differ¬ 
ent  depending  on  whether  the  surfaces  are  rough 
or  smooth.  Figure  2  shows  the  sensitivity  of  heat 
transfer  predictions  to  the  near-wall  spacing  for  a 
representative  rotor  geometry  for  both  smooth  and 
rough  surfaces.  The  results  are  given  in  terms  of 
a  reference  y"*".  The  same  definition  is  used  as  was 
used  by  Boyle  in  reference  13.  The  definition  is: 


tfREF  =  0.17yif2e°  9/s°  1  (10) 


Where  yi  is  the  distance  from  the  surface  of  the 
first  grid  line.  Re  is  the  exit  Reynolds  number  per 
unit  of  length.  y£EF  is  only  a  weak  function  of 
the  surface  distance,  s,  and  s  is  taken  as  the  axial 
chord.  The  same  definition  is  used  for  the  reference 
surface  roughness,  hEEF  except  that  yi  is  replaced 
by  heq. 


*ref  =  0.17/ieqRe°  9/s°  1  (11) 

The  calculations  are  for  fully  turbulent  flow. 
The  smooth  wall  results  shown  in  figure  2a  are  not 
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very  sensitive  to  the  near  wall  spacing.  A  greater 
sensitivity  was  shown  in  reference  13  for  a  similar 
comparison.  The  sensitivity  was  greater  when  the 
flow  was  laminar,  or  close  to  transition. 

Figure  2  shows  a  somewhat  greater  sensitiv¬ 
ity  to  the  near  wall  spacing  for  the  rough  surface 
th an  for  the  smooth  one.  The  reference  values  were 
calculated  assuming  a  smooth  surface,  and  the  cal¬ 
culation  depends  on  the  inverse  of  the  friction  fac¬ 
tor.  Thus,  for  rough  surfaces  the  reference  values 
is  lower  than  one  based  on  the  correct  friction  fac¬ 
tor  by  the  ratio  y/{C f)RovGu/ {C /)smooth*  The 
rough  surface  friction  is  not  known  before  the  anal¬ 
ysis  is  done,  and  equation  11  is  useful  for  determin¬ 
ing  the  grid  spacing.  Consequently,  the  rough  wall 
heat  transfer  shown  in  figure  2b  is  more  sensitive 
to  the  near-wall  spacing  than  the  smooth  wall  data 
shown  in  figure  2a. 

The  primary  purpose  of  this  investigation  was 
to  examine  the  effects  of  surface  roughness  on  heat 
transfer  predictions.  Consequently,  grids  with  close 
near-wall  spacing  were  used.  Calculations  were 
typically  done  using  145  x  54  C-grids.  Approxi¬ 
mately  5000  iterations  were  needed  to  insure  con¬ 
vergence.  The  computations  took  about  700  CPU 
seconds  on  a  Cray-XMP.  It  is  expected  that  further 
optimization  of  the  code  and  calculation  procedure 
would  result  in  faster  convergence. 

Figure  3  shows  the  effect  of  the  different 
crossover  models  for  fully  turbulent  flow.  The  cal¬ 
culations  are  for  a  smooth  wall.  The  difference  in 
the  calculated  heat  transfer  among  the  different 
models  is  not  great.  Of  the  two  blending  mod¬ 
els,  equation  9b  is  closer  to  the  Baldwin- Lomax 
crossover  model.  The  Baldwin- Lomax  model  re¬ 
sults  in  a  heat  transfer  prediction  that  has  greater 
fluctuations  between  adjacent  nodes  than  either  of 
the  other  crossover  models.  This  is  consistent  with 
the  abrupt  nature  of  the  crossover  model.  For  some 
cases  the  solution  stability  or  convergence  might  be 
enhanced  using  a  blended  crossover  model. 

The  equivalent  roughness  height,  heq,  is  a 
function  of  the  actual  roughness  and  the  rough¬ 
ness  density.  A  number  of  correlations  have  been 
proposed  to  obtain  heq.  Sigal  and  Danberg^14)  dis¬ 
cuss  various  correlations  for  obtaining  the  equiv¬ 
alent  height.  The  ratio  of  equivalent  to  actual 
roughness  height,  hcq/h}  is  correlated  as  a  func¬ 
tion  of  the  roughness  density  parameter,  A.  Fig¬ 


ure  4  shows  their  results  with  different  correlations 
for  both  two-  and  three-dimensional  roughness.  A 
is  defined  differently  for  the  different  correlations. 
The  definition  used  by  Sigal  and  Dan  berg  was: 

<**•> 

The  definition  used  by  Dvorak^15)  as  interpreted 
by  Simpson^16)  was: 

Adv  =  (126) 

bf 

The  definition  used  by  Dirling^17^  was: 

Ade  =  a(^)  ”  <I2c) 

Sigal  and  Danberg  reasoned  that  three-dimensional 
roughness  has  a  lower  hcq  because  flow  is  able  to  go 
around  as  well  as  over  a  three-dimensional  obstacle. 
Since  different  definitions  were  used  for  A,  Dvorak’s 
correlation  does  not  necessarily  result  in  a  lower 
value  of  heq. 

In  the  correlation  of  Sigal  and  Danberg,  and  in 
the  correlation  of  Dvorak,  for  the  same  roughness 
shape,  A  is  proportional  to  the  heigh t-to-pitch  ratio 
for  two-dimensional  roughness.  However,  for  three- 
dimensional  roughness  A  is  proportional  to  the 
square  of  this  ratio.  Consequently,  even  for  closely 
spaced  roughness,  only  two-dimensional  roughness 
results  in  A  values  in  the  region  where  the  slope  of 
the  correlation  is  positive. 

Hosni  et  al/8^  conducted  experiments  on 
smooth  and  rough  flat  plates.  The  rough  surface 
test  sections  were  formed  from  machined  aluminum 
plate,  and  had  integral  hemispherical  roughness  el¬ 
ements.  They  were  of  diameter  Z?,  and  were  spaced 
one  of  three  distances,  L,  apart  to  form  a  three- 
dimensional  roughness  array.  This  resulted  in  three 
different  roughness  densities.  Table  I  shows  the 
equivalent  height  ratio  for  different  spacings  using 
the  three  definitions  of  A.  Because  A  was  defined 
differently  by  the  different  investigators,  there  is 
less  variation  in  the  equivalent  height  ratio  than 
might  be  expected  from  figure  4.  The  equivalent 
height  ratio  varied  from  about  0.02  to  somewhat 
greater  than  1.0  as  L/D  varied  from  2  to  10. 
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Figure  5  shows  comparisons  of  predicted  and 
measured  heat  transfer  over  a  flat  plate  for  differ¬ 
ent  equivalent  roughness  heights.  Since  the  Navier- 
Stokes  analysis  used  here  was  developed  for  turbo- 
machinery  applications,  the  flat  plate  results  were 
calculated  as  a  thin  uncambered  airfoil  in  a  low 
solidity  cascade.  The  airfoil  had  an  elliptical  lead¬ 
ing  edge  with  a  10:1  radius  ratio.  This  was  done 
to  avoid  any  separation  in  the  leading  edge  region, 
as  might  occur  with  a  circular  leading  edge.  The 
experimental  data  are  shown  for  clarity  as  indi¬ 
vidual  points,  but  were  derived  from  a  curve  fit 
through  data.  As  the  Reynolds  number  increases 
the  predicted  heat  transfer  becomes  greater  than 
the  experimented  data  for  the  smooth  surface.  The 
predicted  heat  transfer  is  greater  by  nearly  10%. 
In  the  Navier-Stokes  analysis  the  solidity  may  not 
have  been  low  enough.  There  was  some  small  ac¬ 
celeration  of  the  freestream  flow  due  to  boundary 
layer  growth.  This  would  increase  the  heat  transfer 
above  that  for  an  isolated  flat  plate. 

Table  II  compares  the  ratio  of  predicted  to 
measured  heat  transfer  for  the  different  roughness 
spacings  at  a  Reynolds  number  of  400000.  Com¬ 
parisons  are  given  for  the  values  of  heq/ h  from 
the  different  roughness  correlations.  Also  shown 
in  this  table  is  the  ratio  of  predicted  to  measured 
heat  transfer  for  the  smooth  surface.  The  analy¬ 
sis  predicts  the  trends  in  the  experimental  data. 
Generally,  the  best  agreement  is  found  when  the 
lowest  value  of  ht<i/h  for  a  given  spacing  is  used 
in  the  analysis.  The  most  recent  correlation  of 
roughness  data  is  that  of  Sigal  and  Danberg.  This 
is  also  the  correlation  that  agrees  best  with  the 
data.  Sigal  and  Danberg  showed  that  spherical 
segment  data  might  deviate  significantly  from  their 
proposed  correlation.  This  is  the  roughness  geom¬ 
etry  that  most  closely  resembles  the  roughness  el¬ 
ements  used  by  Hosni  et.  al.  A  line  using  the 
extremes  of  the  spherical  segment  data  would  be 
about  a  factor  of  3  less  in  equivalent  height  ratio 
than  the  three-dimensional  correlation  of  Sigal  and 
Danberg.  Consequently,  it  is  not  suprising  that  a 
heat  transfer  analysis  using  the  correlation  of  Sigal 
and  Danberg  would  overpredict  the  heat  transfer 
when  hemispherical  roughness  elements  are  used. 

There  is  little  rough  surface  turbine  blade 
heat  transfer  data  available  that  is  suitable  for 
comparisons  with  analytic  predictions.  Suitable 


data  would  have  well  defined  information  on  the 
roughness  characteristics,  and  not  just  the  rough¬ 
ness  height.  Measurement  techniques  suitable  for 
smooth  blades  may  not  be  suitable  for  rough  sur¬ 
faces.  For  example,  Taylor,  Taylor,  Hosni,  and 
Coleman^18)  showed  that  surface  mounted  sensors 
which  do  not  have  the  same  roughness  character¬ 
istics  as  the  blade  itself  may  indicate  heat  transfer 
rates  different  from  that  of  the  blade.  If  the  surface 
of  the  blade  is  artificially  roughened  using  a  low 
conductivity  material  such  as  sand,  this  material 
may  effectively  move  transition  forward  on  the  sur¬ 
face,  but  also  form  a  thermal  resistance.  These  re¬ 
quirements  suggest  that  experimental  heat  transfer 
data  for  rough  turbine  blades  should  include  base¬ 
line  smooth  surface  data.  The  experimental  data 
of  Tarada(19)  showed  that  for  many  cases  there 
was  no  increase  in  heat  transfer  for  rough  blades 
when  the  flow  is  turbulent.  In  these  data,  in  which 
roughness  was  achieved  by  coating  a  smooth  blade 
with  sand,  the  roughness  only  tripped  the  bound¬ 
ary  layer  close  to  the  leading  edge.  The  early  trip¬ 
ping  of  the  boundary  layer  resulted  in  an  increased 
average  surface  heat  transfer.  However,  for  many 
cases  the  heat  transfer  was  lower  in  the  fully  turbu¬ 
lent  region  for  rough  blades  than  for  smooth  blades. 
The  analysis  does  not  predict  this  behavior.  The 
predicted  effect  of  roughness  for  turbine  blades  is 
similar  to  the  predictions  for  a  flat  plate  shown  in 
figure  5. 

In  contrast  to  the  data  of  Tarada,  the  data 
of  Blair  and  Anderson^10*  show  higher  heat  trans¬ 
fer  for  a  roughened  blade  in  the  fully  turbulent  re¬ 
gion.  While  the  roughness  height  was  given,  no 
information  was  given  to  determine  the  roughness 
density.  Sand  was  used  to  roughen  the  surfaces, 
but  the  tests  were  conducted  in  a  much  different 
thermal  environment.  The  significance  of  an  insu¬ 
lating  layer  may  have  been  less  in  this  environment 
than  in  the  tests  of  Tarada.  The  tested  blade  was 
the  rotor  of  a  large-scale  low-speed  annular  cas¬ 
cade.  The  turbulence  intensity  is  high  in  front  of 
the  blade  due  to  the  presence  of  the  upstream  sta¬ 
tor.  Figure  6  compares  predicted  and  measured 
heat  transfer  for  the  smooth  surface.  Two  predic¬ 
tions  are  shown  in  figure  6a.  The  one  in  best  agree¬ 
ment  with  the  data  was  for  a  turbulence  intensity 
of  10%.  The  augmentation  of  laminar  heat  trans¬ 
fer  was  included  until  the  flow  was  fully  turbulent. 
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The  transition  model,  which  is  very  sensitive  to  the 
turbulence  intensity,  is  the  one  described  in  refer¬ 
ence  13.  The  other  calculation  neglects  the  effect 
of  turbulence  with  respect  to  both  transition  and 
augmentation  of  laminar  heat  transfer.  Compar¬ 
isons  between  the  two  calculations  show  the  effect 
of  freestream  turbulence  to  be  large. 

Rough  surface  predictions  are  shown  in  figure 
6b.  Calculations  are  shown  for  different  values  of 
the  equivalent  height  ratio.  Since  sand  was  used  to 
roughen  the  blade,  an  equivalent  height  ratio  near 
unity  is  expected.  In  fact  this  does  give  reasonably 
good  agreement  with  the  data.  The  pressure  sur¬ 
face  heat  transfer  calculation  is  less  sensitive  than 
the  suction  surface  to  the  assumed  value  of  equiv¬ 
alent  height  ratio. 

Figure  7  shows  heat  transfer  predictions  for 
both  smooth  and  rough  surfaces  of  the  turbine  ro¬ 
tor  blade  tested  by  Liebert^.  Only  a  rough  sur¬ 
face  blade  was  tested.  The  primary  purpose  of  the 
test  was  to  verify  a  heat  flux  measurement  tech¬ 
nique.  Measurements  were  made  on  a  1:1  scale 
SSME  high-pressure  fuel  turbine  rotor  blade.  Be¬ 
cause  of  the  relatively  small  size  only  one  sensor 
was  installed  on  each  blade.  Only  three  blades  were 
instrumented  in  the  test.  To  obtain  a  complete 
mapping  of  the  blade  heat  transfer,  several  blades 
could  be  instrumented.  The  surface  roughness  pro¬ 
files  were  measured.  Both  the  surface  roughness 
height,  h ,  and  the  roughness  density,  A,  were  cal¬ 
culated.  The  equivalent  height  ratio,  Aeq/A,  was 
determined  to  be  between  0.2  and  0.4.  Predictions 
are  shown  for  different  height  ratios  as  well  as  for  a 
smooth  blade.  The  tests  were  conducted  in  a  tur¬ 
bine  tester  rig.  This  rig  matched  the  engine  turbine 
inlet  temperature,  but  operated  at  a  lower  pres¬ 
sure  than  the  actual  engine.  The  lower  pressure 
resulted  in  lower  than  engine  Reynolds  number.  If 
the  roughness  height  were  the  same,  the  value  of 
in  the  engine  would  be  higher  than  in  these 
tests.  These  tests,  like  those  of  Blair  and  Ander¬ 
son,  were  also  in  a  high  turbulence  environment. 
Because  of  the  high  turbulence  intensity  and  high 
Reynolds  number,  fully  turbulent  predictions  are 
shown.  The  data  are  in  reasonably  good  agree¬ 
ment  with  the  prediction  for  an  equivalent  height 
ratio  of  0.1.  The  figure  also  shows  that  had  the 
roughness  been  of  the  same  height,  but  of  a  dif¬ 
ferent  density,  the  heat  transfer  might  have  been 


significantly  higher. 

The  fully  turbulent  predictions  shown  in  fig¬ 
ure  7  used  an  averaging  process  near  the  stagnation 
point.  There  is  symmetry  in  flow  properties  on  ei¬ 
ther  side  of  the  grid  line  that  intercepts  the  stagna¬ 
tion  point.  Consequently,  the  calculated  vorticity 
is  nearly  zero  along  this  line.  This  in  turn  results 
in  a  turbulent  eddy  viscosity,  /**,  that  is  nearly  zero 
along  the  grid  line  which  leads  up  to  the  stagna¬ 
tion  point.  For  grid  lines  adjacent  to  this  line,  the 
vorticity  is  large,  as  is  fit,  In  the  analysis  fit  was 
taken  as  the  average  of  the  values  at  five  grid  lines 
in  the  stagnation  region. 

CONCLUSIONS 

Predictions  for  rough  surface  heat  transfer  ex¬ 
hibit  a  greater  sensitivity  to  the  near  wall  spacing 
than  do  those  for  smooth  surfaces.  This  results  in 
the  need  for  closer  near  wall  grid  spacing  for  the 
rough  surface  predictions. 

The  modified  mixing  length  approach  of  Ce- 
beci  and  Chang  can  be  used  to  predict  the  trends  in 
heat  transfer  rates  with  roughness.  The  disagree¬ 
ment  between  the  analysis  and  the  experimental 
data  could  generally  be  accounted  for  by  the  un¬ 
certainty  in  the  assumptions  regarding  roughness. 
The  primary  uncertainty  is  the  equivalent  height. 

The  results  of  this  investigation  illustrate  the 
importance  of  having  a  well  characterized  surface 
definition  for  experiments  used  to  validate  rough 
turbine  blade  heat  transfer  predictions.  Not  only 
should  the  roughness  height  and  density  be  known, 
but  the  thermal  characteristics  of  the  roughness 
should  also  be  known.  For  the  data  to  be  suit¬ 
able  for  comparison  with  analytic  predictions  it  is 
necessary  to  quantify  any  additional  thermal  resis¬ 
tivity  caused  by  the  roughness. 
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Table  I.  -  Equivalent  height  ratios  for  roughness  geometries  of  Hosni  et  al. 
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Table  II.  -  Ratio  of  predicted  to  measured  heat  transfer  for  a 
flat  plate  at  Rex  =  400000.  Data  of  Hosni  et  al. 
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Fig.  1  Offset  in  mixing  length  as  a  function 
of  roughness  height. 
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Fig.  3  Effect  of  crossover  model  on  blade  heat  transfer 
for  smooth  blade 


ROUGHNESS  DENSITY  PARAMETER,  A 


b)  Rough  surface, h*REF=60. 

Fig.  2  Effect  of  near-wall  spacing  on  calculated 
blade  heat  transfer  distribution. 


Fig.  4  Equivalent  roughness  height  ratio  as  a 
function  of  roughness  density  parameter 
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Fig.  5  Calculated  and  measured  heat  transfer  for 
a  smooth  and  rough  flat  plate 
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Fig.  6  Comparison  of  predicted  and  measured  rotor  blade 
heat  transfer.  Data  of  Blair  and  Anderson. 
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Fig.  7  Comparison  of  predicted  and  measured  rotor  blade 
heat  transfer.  Data  of  Liebert. 


11 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Pubhc  repixting  burden  for  ms  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  me  time  for  reviewing  instructions,  searching  existing  data  sources 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  mis  burden  estimate  or  any  other  aspect  of  mis 
Information  including  suggestions  for  reducing  mis  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports  1215  Jefferson 

Davis  Highway.  Suite  1204,  Arlington,  VA  22202^*302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  {Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

June  1991  Technical  Memorandum 

4.  TITLE  AND  SUBTITLE 

Two-Dimensional  Navier-Stokes  Heat  Transfer  Analysis  for  Rough  Turbine 

Blades 

S.  FUNDING  NUMBERS 

WU-505-62-52 

6.  AUTHOR(S) 

RJ.  Boyle  and  K.C.  Civinskas 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

NASA  Lewis  Research  Center 

Cleveland,  Ohio  44135-3191 
and 

Vehicle  Propulsion  Directorate 

U.S.  Army  Research  Laboratory 

Cleveland,  Ohio  44135-3191 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

E-7755 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

National  Aeronautics  and  Space  Administration 

Washington,  D.C.  20546-0001 
and 

U.S.  Army  Research  Laboratory 

Adelphi,  Maryland  20783-1145 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

NASA  TM-1 06008 

AIAA-91-2129 

11.  SUPPLEMENTARY  NOTES 

Prepared  for  the  27th  Joint  Propulsion  Conference  and  Exhibit  sponsored  by  the  AIAA,  SAE,  and  ASME,  Sacramento,  California,  June  24-27,  1991. 

R.J.  Boyle,  NASA  Lewis  Research  Center;  and  K.C.  Civinskas,  Vehicle  Propulsion  Directorate,  U.S.  Army  Research  Laboratory,  Lewis  Research 

Center,  Cleveland,  Ohio  44135.  Responsible  person,  R.J.  Boyle,  (216)  433-5889. 

12a.  DISTRIBUTION/AVAILABIUTY  STATEMENT 

Unclassified  -  Unlimited 

Subject  Category  34 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

A  quasi-three-dimensional  thin-layer  Navier-Stokes  analysis  was  used  to  predict  heat  transfer  to  rough  surfaces.  Com¬ 
parisons  are  made  between  predicted  and  experimental  heat  transfer  for  turbine  blades  and  flat  plates  of  known  roughness. 
The  effect  of  surface  roughness  on  heat  transfer  was  modeled  using  a  mixing  length  approach.  The  effect  of  near-wall 
grid  spacing  and  convergence  criteria  on  the  accuracy  of  the  heat  transfer  predictions  are  examined.  An  eddy  viscosity 
mixing  length  model  having  an  inner  and  outer  layer  was  used.  A  discussion  of  the  appropriate  model  for  the  crossover 
between  the  inner  and  outer  layers  is  included.  The  analytic  results  are  compared  with  experimental  data  for  both  flat 
plates  and  turbine  blade  geometries.  Comparisons  between  predicted  and  experimental  heat  transfer  showed  that  a 
modeling  roughness  effects  using  a  modified  mixing  length  approach  results  in  good  predictions  of  the  trends  in  heat 
transfer  due  to  roughness. 

14.  SUBJECT  TERMS 

Turbine;  Heat  transfer;  Rough  surface 

IS.  NUMBER  OF  PAGES 

12 

16.  PRICE  CODE 

A03 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

NSN  7540-01  -280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239-18 
298-102 


